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ABSTRACT 



We study the effect of galactic outflows on the statistical properties of the Lycv 
forest and its correlation with galaxies. The winds are modelled as fully ionised spher- 
ical bubbles centered around the haloes in an N-body simulation of a ACDM model. 
The observed flux probability distribution and flux power spectrum limit the volume 
, filling factor of bubbles to be less than 10%. We have compared the mean flux as 

a function of distance from haloes with the Adelberger et al. (ASSP) measurement. 
For a model of bubbles of constant size surrounding the most massive haloes, bubble 
C^l . radii of >1.5 h _1 Mpc are necessary to match the high transmissivity at separations 

< 0.5 h Mpc but the increase of the transmissivity at small scales is more gradual 
than observed. The cosmic variance error due to the finite number of galaxies in the 
Oh" sample increases rapidly with decreasing separation. At separations < 0.5 h _1 Mpc our 

estimate of the cosmic variance error is AF ~ 0.3, 30% higher than that of ASSP. The 
difficulty in matching the rise in the transmissivity at separations smaller than the size 
of the fully ionised bubbles surrounding the haloes is caused by residual absorption of 
neutral hydrogen lying physically outside the bubbles but having a redshift position 
similar to the haloes. The flux level is thus sensitive to the amplitude of the coherent 
velocity shear near halos and to a smaller extent to the amplitude of thermal motions. 
, We find that the velocity shear increases with halo mass in the simulation. A model 

5_i ■ where LBGs are starbursts in small mass haloes matches the observations with smaller 

bubble radii than a model where massive haloes host the LBGs. If we account for the 
uncertainty in the redshift position of haloes, a starburst model with a bubble radius 
of 1 h _1 Mpc and a volume filling factor of 2% is consistent with the ASSP measure- 
ments at the l-1.5cr level. If this model is correct the sharp rise of the transmissivity 
at separations < 0.5 h _1 Mpc in the ASSP sample is due to cosmic variance and is 
expected to become more moderate for a larger sample. 

Key words: cosmology: theory - gravitation - dark matter -baryons- intergalactic 
medium 
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1 INTRODUCTION 

Stellar feedback, in the form of ionising radiation and me- 
chanical energy in galactic winds/outflows, plays an impor- 
tant role in determining the physical properties of galaxies 
and the Intergalactic Medium (IGM). The galactic contribu- 
tion to the hydrogen ionising UV radiation may well exceed 
that of QSOs. Galactic superwinds regulate star formation, 
enrich the IGM with metals, and alter the dynamical and 
observational properties of the IGM close to galaxies. 

Galactic winds are driven by the mechanical energy 
produced by supernova explosions of massive stars. They 



have been studied extensively in nearby starburst galaxies 
(e.g. Heckman et al. 1990). The large velocity shifts between 
stellar and interstellar lines in the spectra of Lyman-break 
galaxies (LBGs) are strong evidence that galactic winds are 
also present in high-redshift galaxies (Pettini et al. 1998, 
2000, Heckman et al. 2000). The inferred wind velocities are 
several hundred kms^ 1 , far above the sound speed of the 
IGM so that strong shocks are likely to be associated with 
these winds. These shocks will heat up, collisionally ionise, 
and sweep up the surrounding IGM near galaxies (Heck- 
man 2000). Adelberger et al. (2003, hereafter ASSP) have 
analyzed a sample of high resolution QSO spectra with lines- 
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of-sight (LOS) passing very close to LBGs. They detected a 
significant decrease of the Lya absorption near galaxies. As 
discussed by ASSP and Kollmeier et al. (2003) the increased 
UV flux close to LBGs falls short of explaining the decreased 
absorption by a factor of a few. The decreased absorption is 
thus taken as evidence for the existence of dilute and highly 
ionised gas bubbles caused by wind-driven shocks. 

Initially a shock propagating in the IGM expands adi- 
abatically. During this phase the gas encountered by the 
shock is heated and collisionally ionised to become part of 
a bubble of dilute and hot plasma. When radiative losses 
become important, a thin cool shell forms at the shock 
front. The shell is driven outwards by the hot low density 
inner plasma (e.g. Weaver et al. 1977, Ostriker & McKee 
1988, Tegmark et al. 1993). Eventually the expansion of the 
plasma comes to a halt when pressure equilibrium with the 
surrounding IGM is reached. To explain the lack of Hi close 
to the haloes, the bubble has to survive for a sufficiently 
long time. We defer the discussion of the conditions for the 
survival of the bubble to another paper. Here we assume 
that the bubbles are long-lived and study the implications 
of long-lived bubbles devoid of neutral hydrogen on QSO 
absorption spectra. 

We will investigate how the correlation between the flux 
in the Lya forest and the distance of the LOS to the closest 
galaxy depends on physical properties of the ionised bub- 
bles around galaxies. For this purpose we use a N-body 
simulation of collisionless dark matter (DM). The basic ob- 
served properties of the Lya forest can be reproduced with 
such simulations using a simple prescription that relates the 
distribution of gas and dark matter in the IGM (Hui & 
Gnedin 1997). With DM simulations it is possible to probe 
a wider dynamical range than with hydrodynamical sim- 
ulations. This allows a better modelling of the statistical 
properties of the IGM-halo relation. 

Measures of the halo-flux correlation, defined for exam- 
ple as the mean flux in the Lya forest as a function of the 
distance to the nearest galaxies, are sensitive to the size of 
the bubbles. Absorption spectra probe the gas distribution 
in redshift space. Thermal motions and coherent peculiar ve- 
locities therefore strongly affect the halo-flux correlation at 
small separations (e.g. Kollmeier et al. 2003). Previous mod- 
elling of the impact of winds on the Lya forest had difficul- 
ties to match the large decrease of the absorption measured 
by ASSP at comoving separation s < 0.5 h _1 Mpc from the 
galaxies (Croft et al. 2002, Bruscoli et al. 2002, Theuns et 
al. 2002, Kollmeier et al. 2003). We will investigate here the 
role of peculiar velocities, thermal motions and the scaling 
of bubble size with halo mass in more detail. We will also 
perform a detailed assessment of the errors to obtain a more 
robust lower limit on the size and thus the volume filling 
factor of the bubbles. 

The paper is organized as follows. In §2 we describe the 
model of the Lya forest which we have used to calculate syn- 
thetic spectra from the DM simulations. In §3, we describe 
the simulation and the corresponding halo catalogs, compute 
three-dimensional flux maps of the simulation, and discuss 
the flux probability distribution function and power spec- 
trum in the presence of winds. The main body of the paper 
are §4 and §5 where we study in detail the impact of galac- 



tic winds on the halo-flux correlation and the corresponding 
observational errors. We conclude with a discussion of the 
results in §6. 



2 THE LYMAN a FOREST AND THE 
MODERATE DENSITY IGM 

The Lya forest is now widely believed to originate from an 
undulating warm ( 10 4 K) photoionized IGM. The Lya opti- 
cal depth in redshift space due to resonant scattering can be 
expressed as a convolution of the real space Hi density along 
the line of sight with a Voigt profile H (Gunn & Peterson 
1965, Bahcall & Salpeter 1965), 



n H1 (x)H [w — x — v p (x),b(x)^ dx . (1) 

-oo 

where <tq — 4.45 x 10~ 18 cm 2 is the effective cross-section 
for resonant line scattering and H(z) is the Hubble constant 
at redshift z, x is the real space coordinate, v p (x) is the 
line of sight component of the Hi peculiar velocity field, Ti. 
is the Voigt profile, and b(x) is the Doppler parameter due 
to thermal/turbulent broadening. For thermal broadening 

6(x) = 13 ( ) ^ kms -1 , where T is the temperature of 
the gas. 

Both x and w are measured in km s~ 1 . The Voigt profile 
is normalized such that f 7i = 1 and can be approximated 
by a Gaussian for moderate optical depths, TL = 1/ (n 1 / 2 ^ x 
exp { — \u) — x — v p (x)] 2 /6 2 } . 

The normalized flux obtained by "continuum fitting" of 
the observed spectrum is related to the optical depth along 
the line of sight as, 

%)^»bsW/U=e- TW , (2) 

where r is the optical depth, w is the redshift space coor- 
dinate along the line of sight, J b s is the observed flux and 
^cont is the flux emitted from the source (quasar) that would 
be observed in the absence of any intervening material. 

On scales larger than a filtering scale which is related 
to the Jeans scale , xj, the IGM traces the DM distribu- 
tion very well. For moderate overdensities the balance be- 
tween adiabatic cooling and photo-heating of the expanding 
IGM establishes power-law relations between temperature, 
neutral hydrogen density and total gas density (Katz et al. 
1996, Hui & Gnedin 1997, Theuns et al. 1998), 

T g oc p| and n HI oc Pg , (3) 

where the parameter /3 is in the range — 0.62, a = 2 — 0.7(3, 
and p g and T g are the density and temperature of the gas, 
respectively. 

On scales smaller than the Jeans scale pressure domi- 
nates over gravity. The gas pressure smoothes the gas dis- 
tribution compared to the distribution of collisionless dark 
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Table 1. The principle parameters of the N-body simulation. 
The number of particles, N, and the particle mass, m, in the 
high resolution region are also listed. 
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0.9 


12067979 9.52 



matter (e.g. Theuns, Schaye & Haehnelt 2000). We thus as- 
sume that p g (x) oc p® m (x), where p° m (x) is the dark mat- 
ter density smoothed in a way that mimics the pressure ef- 
fects. We adopt the smoothing used in SPH simulations (e.g. 
Springel, Yoshida & White 2001) to estimate the gas density 
from the dark matter particle distribution. The Hi density 
n m and the smoothed dark matter field p° m are then related 
as 



n HI (x) 



(4) 



where p| m is the (volume) average of p s dm , and n HI is the 
Hi density at p dm = p B dul . 



3 SIMULATING GALAXIES AND THE FLUX 
DISTRIBUTION 

3.1 The DM simulation 

One of our main goals is to assess the error in the observed 
galaxy flux correlation. This requires a simulation large 
enough to contain a substantial number of DM haloes which 
can be identified as hosts of LBGs. At the same time the sim- 
ulation should still reproduce typical absorption systems. 
This favours the use of a DM simulation for which a larger 
dynamical range can be achieved rather than a computa- 
tionally more expensive hydrodynamical simulation. We use 
a ACDM simulation run with the N-body code GADGET 
(Springel, Yoshida & White 2001). The cosmological param- 
eters are listed in Table 1. It is part of a series of simulations 
of increasingly higher resolution which zoom in on a spher- 
ical region of initial comoving radius 26 h _1 Mpc selected 
from a simulation with box size 479 h _1 Mpc (cf. Stoehr et 
al. 2002). The inner region of the simulation has thus higher 
resolution than the outer regions. We only use the particle 
distribution inside a box of size L = 30 h _1 Mpc which be- 
longs to the spherical region of high resolution. The particle 
mass in this high-resolution region is m = 9.52 x 10 s M /h 
(cf. Table 1). The large volume makes the simulation well 
suited for a direct comparison with the observed galaxy-flux 
correlation. We use a redshift output z = 3 to facilitate a 
comparison with the ASSP measurements. 



3.2 DM haloes as host of LBGs 



is still somewhat uncertain. When LBGs were first discov- 
ered, a simple picture was advocated in which the mass 
of the DM haloes hosting the LBGs scales approximately 
linearly with the luminosity of LBGs, LBGs are long-lived 
and each DM halo hosts one LBG (e.g. Steidel et al. 1996; 
Adelberger et al. 1998, Bagla 1998; Haehnelt, Natarajan 
& Rees 1998). In this "massive-halo" scenario the masses 
of the host haloes are rather large, M J>10 12 M Q /h, and 
the high-redshift LBGs are the progenitors of today's mas- 
sive and luminous galaxies. The rather slow decrease of the 
LBG space density with increasing redshift motivated an 
alternative picture, where most of the observed LBGs are 
interacting starburst galaxies in the process of assembling 
(e.g. Lowenthal et al. 1997, Trager et al. 1997, Kolatt et al. 
1999, Somerville et al. 2001). We will call this scenario the 
"starburst" picture in the following. It predicts that most 
of the observed LBGs are hosted by haloes of smaller mass 
M ~ 10 11 MQ/h, which will eventually merge to form more 
typical galaxies at z ~ 0. Both scenarios predict similar 
clustering properties and cannot be ruled out based on the 
available clustering data (e.g. Wechsler et al. 2001). Direct 
observational constraints on the masses of LBGs are also 
still rather weak. Studies of rest-frame optical nebular emis- 
sion lines imply halo masses in the range 10 10 — 10 11 Mq/1i 
(e.g. Pettini et al. 1998), while near-infrared imaging surveys 
combined with stellar population synthesis models and ex- 
tinction yield virial masses of M ~ 10 11 — 10 12 M /h (e.g. 
Shapley et al. 2001). We will here first adopt the massive 
halo picture i.e. we assume that a dark matter halo con- 
tains only one LBG (e.g. Adelberger et al. 1998) and pick a 
minimum mass M m i n such that the number density of haloes 
M > Mmi n in the simulation is similar to the observed num- 
ber density of high-redshift LBGs. In Section §5.4 we will 
discuss what happens if we relax these assumptions. 

DM haloes were identified with a Friends-of-Friends 
group finding algorithm. Only groups containing at least 20 
dark matter (DM) particles were classified as haloes. Since 
the DM particle mass is 9.52 x 10 s MQ/h, the minimum 
halo mass is about 2 x 10 10 M /h. Fi gure 1 shows a slice 
of thickness Ah = 0.21 h _1 Mpc at z = 3 in the simula- 
tion. The centers of the open circles show the location of 
haloes with M > 5 x 10 10 M Q /h. The radii of the circles are 
proportional to their masses. 

In the simulation the number density of haloes with 
M > 10 11 M Q /h is about 0.04 h 2 Mpc~ 3 , whereas for haloes 
with M > 10 12 M /h it is about 2 x 10~ 3 h 2 Mpc -3 (cf. 
Table 2). We have chosen M min = 5x 10 11 M /h as our ref- 
erence value but have also studied alternative choices. The 
corresponding number density is 4 x 10~ 3 h 2 Mpc~ 3 , which 
is about the number density of LBGs as determined for a 
ACDM cosmology at z = 3 from high-redshift surveys with 
limiting magnitude 25 < 1Z < 27 (Giavalisco & Dickinson 
2001). It should be noted that the halo cumulative func- 
tion fih(mh) as shown in Table 2 is consistent with a Sheth- 
Tormen mass function computed for a ACDM Universe with 
identical parameters as our simulation (cf. Table 1). 



We aim at extracting halo catalogs with statistical prop- 
erties similar to those of observed Lyman-break galaxies 
(LBGs) . However, the relation between LBG and DM haloes 
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Table 2. Number densities of haloes with mass greater than the 
listed four values. 



m h /[10 10 Mq /h] 5 10 50 10 2 

"h(> m h )/[10" 3 h 3 Mpc' 3 ] 81 36 3.9 1.6 
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Figure 1. The particle distribution at z = 3 in a slice of thickness 
Ah = 0.21 h- 1 Mpc. Haloes with masses M > 5 X 10 10 Mq /h arc 
represented by the open circles with radii proportional to their 
masses. 



3.3 Calculating synthetic spectra from the dark 
matter distribution 

3. 3. 1 From density to flux 

The Hi and the flux distributions are calculated from the 
DM distribution which was smoothed with a SPH kernel 
containing 20 particles within the smoothing length. The 
smoothed DM density and velocity fields, p° m and v^ m , were 
interpolated for randomly chosen LOS through the simula- 
tion box. The pixel size of the spectra is A = 60 h _1 kpc 
and the comoving length 30 h _1 Mpc. From the density and 
velocity fields along the LOS we computed the optical depth 
as a function of the redshift space coordinate w = x + v p , 
where v p is the component of v^ m in the line of sight direc- 
tion. Eq. (1) then takes the following discrete form, 

T(w i )=A(z)Aj2[p s d Jxi)] a 'Hii, (5) 

3 

where A = caonm{pd„ 1 )/H(z) and we assume that p B din is 
normalized such that its mean is unity and Hij = H[wi — 
Xj — v p (xj),b(xj)] is a normalized Gaussian profile. 



3.3.2 The temperature distribution 

As we are using a DM simulation we still have to specify a 
temperature. For this we use the fact that temperature and 
density are tightly coupled. At moderate densities a simple 
power- law relation holds, T g ~ p^J 3 (cf. eq 3). In high- 
density regions, line cooling is efficient and the gas rapidly 
cools to T ~ 10 4 K (e.g. Theuns et al. 1998, Dave et al. 
1999). We will thus take the following relation between tem- 
perature and density 

T e =\ f (kY (6) 
10 4 K & > 50 

P 

In the high-density regions where shock heating is impor- 
tant this is only a very crude approximation of the T — p 
diagram of hydrodynamical simulation (Theuns et al. 1998, 
Dave et al. 1999, Croft et al. 2002). However, it is mainly the 
moderate density regions which are responsible for the Lya 
forest. Eq. (6) should thus be sufficient for our purposes. 
The gas temperature T at fixed overdensity is expected to 
evolve with redshift. Schaye et al. 2000 and Theuns et al. 
2002 claim to have detected a peak at z = 3 with decreasing 
temperature towards lower and higher redshifts. Note that 
we have ignored a possible redshift dependence and have 
assumed f = 15000 K ~ 10 4 2 K independent of redshift. 



3.3.3 The effect of wind bubbles 

We have adopted a very simple model for the effect of wind 
bubbles. We simply assume that shocks produce long-lived 
fully ionised spherical bubbles around galaxies (see appendix 
A for a discussion of the physical properties of expanding 
wind bubbles). We assume that the neutral hydrogen Hi 
density is zero in a spherical region of radius r w centered on 
the DM haloes which we have chosen to identify as LBGS. 
Note that this simple model has also been discussed in Croft 
et al. (2002), Kollmeier et al. 2003 (2003a,2003b) and Wein- 
berg et al. (2003). In the simulation, most of the haloes lie 
along filaments, or at the intersection of filaments, and the 
gas distribution around the haloes is often anisotropic on 
the scale of galactic winds (e.g. Croft et al. 2002, Springel & 
Hernquist 2003) . Theuns et al. (2002) have also shown that 
winds propagate more easily into the low density IGM than 
in the filaments, giving rise to highly ionised regions which 
are generally not spherical. It should thus be kept in mind 
that our assumption of a spherical wind bubble will only 
be a reasonable approximation for strong isotropic winds. 
Furthermore the bubbles do not necessarily live longer than 
a Hubble time. A discussion of the physical properties of 
shocks in the IGM can be found in Appendix A3. 



3.3.4 Properties of the synthetic spectra 

Once we have smoothed the dark matter density p s dln , the 
flux distribution depends only on A and a and the size dis- 
tribution of wind bubbles. We have assumed a = 1.6 and 
constrain A from a random sample of 500 lines of sight of 
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Figure 2. The flux power spectrum (left panel) and PDF (right panel) at z = 3 with wind bubbles of radius r w = (solid line), r w = 0.5 
(long-dashed line), r w = 1.0 (short-dashed line) and r w = 1.5 h _1 Mpc (dotted line). The results are compared to the observations of 
McDonald et al. (2000), shown as filled symbols. The triangles arc the various volume averaged statistics in real space. The small plot 
in the right panel shows the PDF for filling factors of zero and ten percent (cf. text). The wavenumber k is in unit of kms -1 . 



comoving length L — 30 h 1 Mpc by demanding that the 
mean flux (F) is equal to the observed value {F) — 0.67 
(Rauch et al. 1997, McDonald et al. 2000). The effect of 
varying the assumed size distribution of wind bubbles are 
discussed in the following sections. The spectral resolution 
is At> ~ 6 kms -1 , somewhat smaller than in the ASSP data, 
where it is Aw <^10 kms~ . Note that the mean inter-particle 
distance in the simulation is ~ 0.2 h _1 Mpc. This is sufficient 
to resolve F on a scale r ~ 0.25 h _1 Mpc ~ 27 kms -1 as m 
the observations, but is somewhat too low to resolve all the 
features of the Lya absorption lines (e.g. Theuns et al. 1998, 
Bryan et al. 1999). We also add noise with signal-to-noise 
S/N=50 per pixel of width A7r = 2.5 kms~\ 



3.4 Volume averaged statistics: The flux PDF and 
the flux power spectrum 

We will first investigate some basic flux statistics of our 
synthetic spectra to check if they are a reasonable repre- 
sentation of observed spectra and to study the effect of the 
wind bubbles. We begin with the flux probability distribu- 
tion function (PDF) defined as the fraction of the volume 
of the IGM which has a flux value within a certain range. 
The right panel of Fig. 2 compares the flux PDF at z — 3 
for wind bubbles with different radii r w . The filled circles 
show the observed PDF from McDonald et al. 2000, while 
the solid, long-dashed, short-dashed and dotted curves cor- 
respond to r w = (i.e., no wind), 0.5, 1 and 1.5 h _1 Mpc, 
respectively. Following McDonald et al. (2000) we use 21 
bins of width AF = 0.05 to compute the flux PDF. The 
PDF of spectra without wind bubbles (solid curve) agrees 
very well with the PDF of the observed spectra (filled cir- 
cles). Including wind bubbles with r w = 0.5 h _1 Mpc has no 



significant effect on the flux PDF, and even when r w is as 
large as 1.5 h~ Mpc there is only a slight discrepancy with 
the flux PDF of the observed spectra in saturated regions 
with F ~ 0. The small effect of the wind bubbles is due to 
their small volume filling factor /hii, which is about 5% for 
the model with r w = 1.5 h~ Mpc. This is consistent with the 
findings of Theuns et al. (2002) and Weinberg et al. (2003). 

We can place an upper bound on the volume filling fac- 
tor, /hii, of bubbles by demanding a good match between 
the simulated and observed PDFs. We have computed the 
flux PDF for several bubbles radii for which /hii is in the 
range 0.01-0.5. To account for the uncertainty in the mass of 
LBGs, we performed three computations for each /hii : wc 
selected DM haloes with mass M larger than 5 x 10 10 , 10 11 
and 5 x 10 11 M /h, and adjusted r w which gives the desired 
/hii- We find that independent of halo mass, a filling factor 
/hii <;0.1 yields a poor match to the observations, especially 
in the range 0.2 <J0.8, as shown in the small plot inside 
the right panel of Fig. 2. 

In the left panel of Fig. 2 we show the dimensionless 
one-dimensional flux power A|(fc) as a function of k (in 
kms -1 ), defined as 

A 2 F (k) = -k\5 k \ 2 , (7) 

where Sk is the Fourier Transform of the flux contrast 
(F/{F) — 1). The flux power spectrum was calculated with 
a FFT routine from 500 synthetic spectra. Note that our 
simulation box and therefore the spectra are not periodic 
which may lead to small artifacts at small scales. The flux 
power of our synthetic spectra and that of the observed spec- 
tra obtained by McDonald et al. (2000) also agree very well. 
Bubbles of size as large as r w = 1 h _1 Mpc have again a very 
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a/ h _1 Mpc s / h _1 Mpc s / h _1 Mpc 



Figure 3. Left panel : Contours of the 2D conditional flux F(tt, <t) around the nearest LBG-size halo in the simulation at z = 3 (cf. text). 
The (comoving) bin size is As = 0.25 h — J Mpc. The shaded pixels have a flux within 0.05 of the mean flux level, (F) = 0.67. Middle 
panel : The mean transmitted flux F(s) at distance s = Vvr 2 + a' 2 , compared to F(0, a). The dotted curve shows the mean flux level. 
Right panel : The mean transmitted flux, when LBG-size haloes were surrounded by wind bubbles with radius r w = 0.5, r w = 1 h _1 Mpc 
and r w = 1.5 h _1 Mpc, respectively, as indicated in the figure. The Adelberger et al. (2003) measurements are shown as filled symbols. 
The shaded area is the lex error for the model r w = 0.5 h _1 Mpc, and was computed from N(< s) shown as an histogram (right axis). 



small effect due to their small volume filling factor (e.g. Croft 
et al. 2002, Weinberg et al. 2003). 

To demonstrate the role of peculiar velocities, the tri- 
angles in Fig. 2 show the flux PDF and the power spectrum 
in real space for r w = obtained by setting the peculiar 
velocities and thermal motions to zero. The real space flux 
PDF and power spectrum significantly deviate from the ob- 
servations. 



4 PROBING CORRELATIONS BETWEEN 
FLUX AND HALOES /LBGS 

4.1 Statistical measures of the flux-halo 
correlation 

There is a variety of statistical measures which can describe 
the correlation between the flux along a LOS and nearby 
galaxies. Here we will focus on the conditional flux function, 
F. F is the mean flux at real-space distance |s| = s to the 
next halo, 

F(s) = ^^F(s|halo), (8) 

where Nh is the total number of haloes, F(s|halo) is the 
value of the flux at a distance s from the nearest halo. Note 
that we use the position of the haloes which we identify 
as host haloes of LBGs as a proxy for the galaxy position. 
Results will thus depend on how this identification is done 
and we will later explore several possibilities. The summa- 
tion is performed over a representative sample of haloes. The 
function F can be expressed in terms of the unconditional 
two-point correlation £ h f as 



F(s) = (F>[1 + Cm(s)] . (9) 

Relation (9) also holds in redshift space. In the following tv 
and a are the coordinates of s parallel and perpendicular to 
the LOS, respectively. 



4.2 Calculating the mean flux as function of the 
distance to the nearest halo 

For the numerical simulation we know the full 3-dimensional 
flux field. We could therefore estimate F(s) from the trans- 
verse correlation F(tt — 0, a). However, real observation will 
sample the flux field F(s) more sparsely. ASSP have there- 
fore estimated F(s) from the averaged 2D halo flux corre- 
lation function £hf using eq. (9) in order to maximize the 
signal. We proceed as ASSP do and calculate F(s) as fol- 
lows. We first pick a LOS direction. We then produce an 
ensemble of spectra along this direction with random off- 
sets perpendicular to the chosen direction. For each pair of 
halo/spectrum we determine the flux F(tt, a), the distance a 
at the point of closest approach, and the distance n between 
the halo and a given pixel of the spectrum. This process is re- 
peated for an ensemble of randomly chosen LOS directions. 
We then obtain F(s) by annular averaging F(tt, a), and bin- 
ning in bins of comoving size As. We further average over all 
haloes and LOS. In the left panel of Fig. 3 we plot the full 
2D distribution of the average flux F(n,a) around haloes as- 
sumed to host LBGs. The haloes were assumed to have no 
galactic wind bubbles. The bin size is 0.25 h _1 Mpc, twice 
smaller than that we use for F. The shaded area shows bins 
for which F is within 0.05 from the mean flux (F). Contours 
are for flux levels increasing from 0.2 to 0.6 with decreasing 
line width. The contours are compressed along the line of 
sight as a result of the peculiar and thermal motions of the 




Figure 4. Left and Middle panels : The transmitted flux in rcdshift space (solid), with gas and halo velocities set to zero but with 
non-zero Doppler parameter b (long-dashed), and with velocities and b set to zero (short-dashed). Panel to the right : The transmitted 
flux as a function of the mean IGM temperature T as indicated on the figure. In each panel the horizontal dotted line shows the mean 
flux level, (F) = 0.67. 



gas. At small separation the transmissivity is much smaller 
than the mean due to the increased density near haloes. 

Redshift distortions cause F(s) (as computed in ASSP) 
to be different from the transverse correlation F(n — 0,a). 
This can be seen in the middle panel of Fig. 3. The effect is 
important for s ~ 1 — 3 h _1 Mpc. Note that previous work 
compared also the ASSP measurements with F(s) obtained 
from numerical simulations (e.g. Croft et al. 2002, Kollmeier 
et al. 2003, Bruscoli et al. 2003). 



4.3 Estimating errors due to sparse sampling 

The sample of ASSP contains 431 galaxies in 7 fields con- 
taining a QSO (one field contains two QSOs). At large scales 
(> a few Mpc) the number of galaxy LOS pairs is reason- 
ably large. However at small scales the number of available 
galaxies becomes very small. There are e.g. only 3 galaxies 
within a transverse separation of 0.5 h _1 Mpc to the line-of- 
sight to a QSO. Obviously Poisson errors will then be large 
and estimating these errors is thus important. In the follow- 
ing we use the term cosmic variance to describe the error 
resulting from the finite number of galaxies and QSO spec- 
tra. At small scales cosmic variance is the dominant source 
of error. Other sources of errors also affect the calculation of 
F, mainly the uncertainty of the mean flux {F}, and the red- 
shift position of the galaxies. These errors will be discussed 
in Section §5.3. 

About one LOS from a bright QSO is expected to pass 
through a region of angular dimensions comparable to that 
of our simulation box at z = 3. However, the correlation 
length of the flux is around a few h _1 Mpc. There will thus 
be about a hundred LOS through the simulation box which 
represent independent flux measurements. We therefore es- 
timate the cosmic variance error as follows from the simu- 



lation. For each separation bin of width As = 0.5 h _1 Mpc 
we perform Monte-Carlo realizations of F(s), where F(s) 
is averaged over the halo/spectrum pairs lying in the given 
bin. The cosmic variance error is then the la scatter around 
F(s). The error of F(s) depends on the total number of 
galaxies/haloes with separation s. A sample of about 40 LOS 
through our simulation box gives a cumulative halo function 
JV(< s) close to the cumulative galaxy function Albg of the 
ASSP sample for separation s < 1 h _1 Mpc. Note also that 
the number of haloes with M > M m i n in the simulation is 
only about 110. This means that we severely oversample the 
halo catalog when we estimate the cosmic variance on large 
scales. However, at the small scales in which we are mainly 
interested this should have little effect. 



5 THE HALO FLUX CORRELATION OF 
OBSERVED AND SIMULATED SPECTRA 

5.1 Wind bubbles and the halo flux correlation 

We first explore the shape of F assuming that all wind bub- 
bles have the same radius, r w , independent of the halo mass 
that they surround. In the right panel of Fig. 3 we plot F for 
haloes with mass M > M m i n . We show curves for r w = 0.5, 1 
and 1.5 h _1 Mpc as indicated in the figure. The correspond- 
ing filling factors are /hii ~ 0.2%, ~ 1.6% and ~ 5.3% 
respectively, small enough to ensure that the flux PDF, the 
power spectrum and the line distribution of the Lya forest 
are not affected. The dashed-dotted curve corresponds to F 
in the absence of outflows. The filled symbols show the ob- 
served correlation as given in ASSP. The shaded area shows 
the cosmic variance error, and is calculated (cf. Section §4.3) 
from the cumulative halo distribution JV(< s) also shown as 
a histogram in the right panel of Fig. 3 (right axis). The 
cosmic variance error is only shown for r w = 0.5 h _1 Mpc. 
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Figure 5. left panel : The velocity dispersion along the line of sight a v near haloes of mass M > 10 11 Mg/h. The solid line is a v averaged 
over bins of size AlogM = 0.25. Middle and right panels : the transmitted flux at separation s < 0.5 h -1 Mpc, Fq, as a function of M. 
The horizontal dashed line is (F) = 0.67. 
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Figure 6. One-dimensional gas density, velocity (real space) and transmitted flux (rcdshift space) along two lines of sight extracted 
from simulation at z = 3, running trough the centers of haloes with mass M = 2 X 10 11 (left panel) and M = 3 X 10 12 M /h (right 
panel). The halo positions along the line of sight are shown by vertical lines. The transmitted flux is computed with winds (solid line) 
and without winds (dotted line). The bubble radius is r w = 0.5 h _1 Mpc. 



The trend of decreasing mean flux with decreasing dis- 
tance to the next galaxy in the observations is well repro- 
duced by the numerical simulations. At small separation 
this trend is inverted if the wind bubbles are included. As 
Kollmeier et al. (2003a, 2003b) we find that the strong in- 
crease of the mean flux level at s :>0.5 h _1 Mpc, is only re- 
produced for bubbles with size as large as r w J>1.5 h _1 Mpc. 
We thus confirm the difficulty of other authors to reproduce 
the mean flux at small separation unless rather large bubbles 
are included. Bubbles of this size are, however, difficult to 



explain even with very efficient galactic superwinds (Croft 
et al. 2002). 

There are thus two possibilities: the simple model of 
spherical wind bubbles of constant size around massive 
haloes does not describe the effect of winds of observed 
galaxies properly or the observed strong increase of the flux 
level in the ASSP is a statistical fluke. We will explore both 
possibilities in the following sections in more detail. 
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Figure 7. Left panel : The transmitted flux F for bubbles of size r w = 0.5 h Mpc after adding a Gaussian deviate of variance cr z = 0, 
a z = 150 and 300 kms -1 to the LOS halo position. Middle panel : The cumulative halo function (histogram) and the la error (shaded 
area) for the model with a z = 150 kms -1 . Right panel : The transmitted flux with a mean flux level (F)=0.62, 0.67 and 0.72. The filled 
symbols are the measurements by Adelberger et al. . 



5.2 The effect of peculiar motions and thermal 
broadening 

Kollmeier et al. (2002) have pointed out that the peculiar 
velocities of matter infalling onto the galaxy/DM halo ex- 
pected from hierarchical structure formation should have a 
strong influence on F(s) at small separations. In Fig. 4 we in- 
vestigate the effect of coherent small-scale peculiar velocities 
and thermal velocities. In the left and middle panel we show 
F(s) with and without peculiar and thermal velocities for 
wind bubbles of size r w — 0.5 h _1 Mpc and r w = 1 h _1 Mpc, 
respectively. Without peculiar velocity the increase of the 
mean flux at small separation is much more pronounced. 
This is easy to understand. In redshift space the velocity 
shear of the infalling surrounding material fills in the cav- 
ity which is caused by the fully ionised wind bubble in real 
space. For the smaller bubble size such an effect is also vis- 
ible for the thermal velocities alone. Neutral hydrogen Hi 
lying at the boundary of the ionised bubbles leads to signif- 
icant Lya absorption at the redshift halo position when the 
ratio 6/r w is large enough. The peculiar velocity will, how- 
ever, always dominate over the effect of thermal velocities. 
This can be seen in the right panel where we have scaled 
the temperature of the gas up and down. This has very lit- 
tle effect on F(s) calculated with peculiar velocities even for 
the smaller bubble size. We have also computed F(s) for 
various values of the density threshold in eq. (6). We found 
only small differences for 10 <^8p/p <^100. Such overdensi- 
ties occur typically within <J0.2 h _1 Mpc (comoving) from 
the halo (Kollmeier et al. 2003) which is smaller than the 
bubble radii considered here. Note that the good agreement 
of the flux PDF of the synthetic spectra with the McDonald 
et al. (2000) data found in section 3 disappeared if peculiar 
velocities were set to zero. This indicates that our simula- 
tions have similar peculiar velocities as the gas responsible 
for the observed absorption - at least in a volume-averaged 
sense. 



In the left panel of Fig. 5 the black dots show the 
velocity dispersion for haloes of mass M. For each DM 
halo the velocity shear a v along the line of sight was cal- 
culated along three perpendicular LOS directions for an 
interval An = 2 h~ Mpc centered on the halo position. 
The histogram shows the mean velocity dispersion aver- 
aged over bins with A log M = 0.25. The velocity shear 
clearly increases with halo mass. This trend is present for 
a range of An, 1 h _1 Mpc <.An <^2 h _1 Mpc. The mid- 
dle and right panels of Fig. 5 show the transmitted flux 
at separation s < 0.5 h~ x Mpc, Fo = F(s < 0.5), as a 
function of the halo mass M assuming that all haloes with 
M > 10 11 Mq/1i are surrounded with bubbles of radius 
r w = 0.5 and r w = 1 h _1 Mpc, respectively. As expected Fo 
is significantly larger for r w = 1 h _1 Mpc. The dependence 
on halo mass appears to be stronger for r w = 1 h _1 Mpc 
than for r w = 0.5 h~ Mpc. In the former case, Fo decreases 
with increasing M by about AF = 0.20 from M = 10 11 
to M = 10 12 M /h. This is consistent with the trend in 
the left panel which predicts that low mass haloes satisfy 
the condition r w ;>(1 + z)a v / H (z) more often than massive 
haloes as a result of the a v — M dependence. For smaller bub- 
bles, r w <J0.5 h _1 Mpc, this appears to be rarely the case, 
independent of the mass of the halo. 

To investigate the effect of the velocity field near haloes 
further we show the real space density and velocity field of 
the gas, as well as the absorption spectrum along two lines 
of sight passing haloes with M = 2 x 10 11 (left panel) and 
3 x 10 12 Mq /h (right panel) in Fig. 6. Solid and dotted line 
are for the case with and without wind bubbles, respectively. 
The radius of the bubbles is r w = 0.5 h _1 Mpc for both 
haloes. The difference between the two cases is striking . 
For the smaller mass halo in the left panel the flux level at 
the position of the halo is F ~ 0.9 if the bubble is included 
whereas F saturates at the redshift space position of the 
more massive halo in the right panel. Since the density fields 
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near the haloes in both panels show similar features, the 
difference in the flux values is mainly due to the difference 
in the velocity fields near the two haloes. The velocity shear 
in the less massive halo (left panel) is quite modest. On 
the contrary, the velocity field near the massive halo (right 
panel) has a large infall on a physical scale of a few hundred 
kms -1 , which moves some of the moderate density peaks 
close to the redshift space position of the halo. This effect 
occurs for a large fraction of haloes, and thereby reduces the 
overall expected F at small scales by a few tens of percent. 



5.3 The role of galaxy redshift error, LOS 
smoothing, and flux normalization 

There is a number of other uncertainties which have to be 
taken into account when comparing the observed F(s) with 
that calculated from the numerical simulation. Most im- 
portant are the uncertain redshifts of the LBGs. There are 
large systematic shifts of several hundred kms -1 between 
the redshift of nebular emission lines, Lya emission and in- 
terstellar absorption lines which make the assignment of a 
redshift somewhat ambiguous. ASSP give Az > 0.002 or 
about f50 kms" 1 as a typical error. We know however rel- 
atively little about the dynamical state of the responsible 
gas and there may well be larger systematic errors. The left 
panel of Fig. 7 shows the effect of adding a Gaussian dis- 
tributed error to the redshift of the DM haloes. The solid 
curve shows the same model as the solid curve in the right 
panel of Fig. 3 ("massive" haloes, r w = 0.5 h _1 Mpc). Long- 
dashed and short-dashed curves are for redshift errors of 
o~ z = f50 kms" 1 and 300 kms" 1 , respectively. As expected 
introducing redshift errors smoothes out the large depres- 
sion of the flux level at distances of up to ~ 2 h" 1 Mpc to 
the line-of-sight. This substantially improves the agreement 
with the ASSP measurements on these scales. 

In the middle panel of Fig. 7 the shaded area shows the 
la error of the simulated F(s) for our ensemble of LBG- 
size haloes with mass M > 5 x 10 M©/h. The error- 
bars show the errors quoted by ASSP. To make the com- 
parison easier the errors are plotted for the model with 
a z — 150 kms" 1 . The solid histogram shows the cumulative 
halo function N(< s) (right axis) as a function of separation. 
At s < 1 h _1 Mpc the cumulative halo function N(< s) is 
close to the cumulative galaxy function iVi.BG of the ASSP 
sample, with N(s < 0.5) = 3 and N(s < 1) = 11 (resp. 3 
and 12 in the ASSP sample) . Our errors are generally larger 
than that of ASSP, at s < 0.5 h" Mpc by about 30 percent. 
At larger scales the difference is smaller but note that there 
our sample is about 20% larger than that of ASSP. Note also 
that the errors decrease faster than those of the transverse 
correlation F(0,a). 

The mean flux level of QSO spectra varies significantly 
from spectrum to spectrum (10-20 percent). We thus also 
investigated how F depends on the value of the mean flux 
{F) . In the right panel of Fig. 7 we show the effect of varying 
(F). The long-dashed curve is our fiducial "massive" haloes 
model, which has r w = 0.5 h~ Mpc and a z = 150 kms" 1 . 
The upper (solid) and lower (short-dashed) curves show 
(F) — 0.72 and 0.62 respectively. A larger (lower) mean 
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Figure 8. Panels to the left : F(s) for a mass-dependent r w 
as indicated in the figure. Panels on the top row have r w = 
0.5 h" 1 Mpc, whereas panels on the bottom row have r w = 
1 h" 1 Mpc. The horizontal dotted line shows the mean flux, 
{F) = 0.67. Panels to the right : the mean transmitted flux F(s) 
for three values of the halo mass limit M m i n . The shaded area is 
the ltr error, attached to the model r w =const. 

transmitted flux (F) raises (lowers) F(s) almost linearly by 
the same amount AF — ±0.05 at all separations. (F) = 0.67 
matches the large scale ASSP measurements best. 



5.4 The effect of varying the bubble model 

5.4.1 F with halo-mass dependent bubble size 

So far we have considered bubbles with sizes that are inde- 
pendent of mass. Unfortunately, it is far from obvious what 
scaling with mass is appropriate. Obviously, on average more 
stars and a larger total energy input from stellar feedback 
are expected in the more massive haloes. However, simple 
binding energy arguments show that winds should escape 
more easily from low mass haloes (Larson 1974, Dekel & Silk 
1986). This may be properly taken into account by a con- 
stant bubble radius but this is very uncertain. We therefore 
investigate the effect of different power-law scalings of the 
bubble size with mass. The left panels of Fig. 8 compare the 
case of constant bubble size with those assuming a power- 
law relation r w (M) oc M v with v = +1/2 and -1/2, at the 
minimal halo mass M m i n = 5 x 10 11 Mq/1i. The top panels 
are for r w = 0.5 and the bottom panels are for r w = 1.0. 
For v — +1/2 F(s) is very sensitive to the value of r w and 
the mean flux at s < 1 h" x Mpc rises to values as large as 
F — 0.8 for the larger of the two bubble sizes improving the 
agreement with the results of ASSP. However, this is at the 
cost of the good agreement with the data on larger scales. 
It appears thus difficult to reproduce the observations of 
ASSP at small and large scales simultaneously by changing 
the mass-radius relation of the bubbles. 
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Figure 9. Left panel : Comparison between the "massive- halo" picture (dashed curve) and the "starburst" picture (solid curve). The 
shaded area shows the la error. A Gaussian deviate of variance a v = 150 kms -1 was added to the halo positions Right panel : A 
comparison between the (real space) two-point correlation function of LBGs haloes in the "starburst" and "massive-halo" models. The 
filled and empty circles show £lbg m the "starburst" and "massive-halo" models respectively. On comoving scale r ~ 0.5 — 5 h — 1 Mpc, 
£lbg m the "starburst" model follows a power-law of index 7 = 1.6 and correlation length ro = 2.5 h -1 Mpc. 



In the right panels of Fig. 8 we investigate the effect of 
changing the minimum mass above which DM haloes are as- 
sumed to be surrounded by wind bubbles. Changing M m i n 
will change the total space density of the halo as given in Ta- 
ble 2. There is less absorption if M m i n decreases, mainly as 
a result of the increasing bubble filling factor . With a bub- 
ble radius r w = 1 h _1 Mpc, the filling factor is fun ~ 15% 
and - 1% for M > 10 11 and M > 10 12 M Q /h, respectively. 



5.4-2 The starburst model 

So far we have assumed that LBGs are hosted by the most 
massive haloes available and have introduced a minimum 
mass to match their space density with that observed for 
LBGs. This assumes that LBGS are long-lived and have the 
same brightness over a fair fraction of the Hubble time. This 
picture may be in conflict with the lack of the evolution of 
the space density of bright LBGs towards very high redshift 
which is more easily explained if LBGs are starbursts with a 
duty cycle of high star formation rate of less than 10%. This 
would mean that there is an about a factor ten or more larger 
number of DM haloes needed to host these LBGs than in the 
massive halo picture and the DM haloes hosting LBGs are 
less massive on average. As should be apparent from the pre- 
vious section, for a given bubble size surrounding a typical 
LBG halo, this should reduce the effect of the velocity shear. 
It is then easier to reproduce the high mean flux levels at 
small separation. We investigate this alternative "starburst" 
picture here in more detail. Kolatt et al. (1999) developed a 
simple, ad hoc procedure to assign LBGs to colliding haloes 
identified in their N-body simulation and were able to re- 
produce both the number density and clustering properties 
of LBGs. Here, we will do something simpler and randomly 



select haloes with mass M > 10 11 Mq/1i such that their 
number density matches the observed LBGs number density, 
n LBG ~ 0.004 h 2 Mpc -3 . This corresponds to a duty cycle of 
ten percent. Note that only 10% of the haloes with a mass 
above 10 11 M /h are haloes as massive as our "massive- 
halo" picture i.e. have masses M > 10 12 M0/I1 (cf. Ta- 
ble 2) . Most of the LBG-size haloes have M <^5x 10 11 M /h 
in the "starburst" scenario. We have produced 120 differ- 
ent Monte-Carlo realizations of such catalogs of starburst- 
ing haloes. In the left panel of Fig. 9 we compare F(s) for 
the "massive-halo" model (solid curve) and the "starburst" 
picture (dashed curve). For both we assume a bubble radius 
r w = 1 h _1 Mpc and a redshift error of the halo position 
of er z = 150 kms -1 . At s <1 h _1 Mpc F(s) is larger in the 
"starburst" model by about 20%. We also plot the la error 
for the starburst model and the ASSP measurements. Our 
starburst model with r w = 1 h _1 Mpc and a volume filling 
factor of 2% appears to be fully consistent with the ASSP 
measurement. Even at s <0.5 h _1 Mpc, where our estimate 
of the cosmic variance error is 30% larger than that of ASSP, 
F(s) in the starburst model is only slightly more than la be- 
low the ASSP measurement despite a difference in the mean 
flux of AF(s) = 0.3. Note that for the quoted volume filling 
factor we did only take into account the bubbles around the 
haloes that host LBGs at a given time. If the bubbles last 
longer than the LBG phase the volume filling factor would 
be higher by the same factor. 



5.4.3 The clustering of LBG haloes 

As discussed above our procedure to identify the host haloes 
of LBGs in the starburst model is rather crude. Unlike Ko- 
latt et al. (1999) and Wechsler et al. (2001) we do not iden- 
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tify "colliding" haloes but choose a fraction of less massive 
haloes at random. It is thus not obvious that the host haloes 
chosen in this way will reproduce the clustering properties of 
LBGs. In the right panel of Fig. 9 we compare the two-point 
correlation function of our "starburst" and "massive-halo" 
models. The correlation function in the starburst model was 
computed from 120 different Monte-Carlo realizations of cat- 
alogs of starbursting haloes. The correlation lengths for the 
"starburst" and "massive-halo" models are 2.5 h _1 Mpc, and 
3.5 h _1 Mpc, respectively. The observed correlation length 
is r ^3-4 h _1 Mpc (Adelberger et al. 1998, Adelberger 
2000, Arnouts et al. 1999, Giavalisco & Dickinson 2001). The 
massive halo picture is thus in somewhat better agreement. 
However, the difference is small and as discussed above, Ko- 
latt et al. (1999) and Wechsler et al. (2001) find that they 
can reproduce the observed clustering properties if they se- 
lect "colliding haloes" as host galaxies of starbursting LBGs. 
We thus do not investigate this issue any further here. 



5.5 Cosmic Variance 

To understand better why the cosmic variance error becomes 
so large at small separation we plot the probability distribu- 
tion P(Fo) in Fig- 10. The distribution is very broad. About 
40% of the haloes have Fq larger than 0.8. The probability 
of finding three (uncorrelated) haloes with F > 0.8 at sepa- 
ration s < 0.5 h _1 Mpc is thus still ~5-10%. The symbols at 
the bottom of Figure 5 show the increase of Fo if the ther- 
mal motions and the velocity shear are artificially reduced 
by a factor of two. The broad probability distribution also 
explains why the suggestion of Croft et. (2002) that LBGs 
are preferentially located in haloes with a low-density envi- 
ronment has such a strong effect on the mean transmission. 



6 DISCUSSION 

We used mock spectra obtained from numerical simulations 
to investigate the effect of fully ionised bubbles around the 
galactic haloes hosting LBGS on the Hi distribution in the 
IGM as probed by the Lyev forest in QSO spectra. By 
matching the probability distribution and the power spec- 
trum of the flux of mock spectra obtained from a N-body 
simulation of a ACDM model with those of observed spectra 
we derive an upper limit of 10% on the volume filling factor 
of bubbles. We further considered the mean transmitted flux 
at a given distance from galactic haloes. To model the ob- 
servational data we computed the mean transmitted flux F 
from an annular average of the two-dimensional conditional 
flux F(-K,a), in the same way as done by Adelberger et al. 
(2003). As previous investigations we find that the observed 
decrease of the transmission at separation 1 < s < 5 h _1 Mpc 
to a halo is well reproduced by the simulated spectra. It 
is due to the increased matter density around haloes. We 
further find that an increase of the transmissivity at the 
smallest separations s < 0.5 h _1 Mpc requires bubble radii 
of r w = 1 — 1.5 h _1 Mpc. This is two to three times larger 
than the separation at which the increase of the transmis- 
sivity is observed. As discussed by Kollmeier et al. 2003 
this is because in redshift space the velocity shear of the 
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Figure 10. PDF of the transmitted flux at separation s < 
0.5 h _1 Mpc, Fo, in the "starburst" picture. The filled circle is 
Fo for s < 0.5 h -1 Mpc as quoted by Adelberger et al. , whereas 
the filled square is the "starburst" model. Empty symbols show 
Fo when the velocity shear and /or the Doppler parameter around 
haloes arc reduced by a factor of two. 



material infalling onto the halo fills in the cavity which is 
caused by the fully ionised wind bubble in real space. For 
a model where the haloes are surrounded with bubbles of 
size r w = 1 — 1.5 h _1 Mpc there is no effect on either the 
flux PDF or the flux power spectrum. The flux PDF and 
the flux power spectrum of our simulated spectra are fully 
consistent with those of observed spectra. This is because 
the volume filling factor of the ionised bubbles is still small. 
We have investigated a variety of scalings of bubble size 
with halo mass. We also compared a model where LBGs 
are long-lived and are hosted by the most massive haloes 
with a model where LBGs are short-lived starburts and are 
hosted by more numerous less massive haloes. Unlike Kolatt 
et al. (1999) and Wechsler et al. (2001), we do not identify 
starbursts as "colliding haloes" but choose a fraction of less 
massive haloes at random. Matching the observed transmis- 
sion at scales 1 < s < 5 h _1 Mpc is not difficult and there is 
certainly more than one way to do so. Introducing redshift 
errors at the expected level improves the agreement with the 
observations. The increase in transmission at small separa- 
tion in our model spectra spectra is, however, generally more 
gradual than that of the observed measurements. As other 
authors we found it challenging to match the observations 
at the smallest separation. Numerical simulations including 
galactic winds and simple analytical estimates as the one in 
Appendix A3 agree that the required bubble sizes are diffi- 
cult to achieve (Croft et al. 2002, Bruscoli et al. 2002). We 
have thus particularly looked into models which are consis- 
tent with the observed measurements and are most "eco- 
nomical" in terms of bubble size. The velocity shear around 
haloes increases with increasing mass. Haloes of lower mass 
require thus smaller bubble radii to explain the observed 
increase of the transmission at s < 0.5 h _1 Mpc. This may 
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favour a model where LBGs are starburst galaxies hosted 
by more numerous less massive haloes. We have taken spe- 
cial care in estimating the errors on the flux near galaxies 
in a similar way as done in the observations. For small sep- 
arations we found that the expected error due to cosmic 
variance is about 30% larger than estimated by ASSP. Our 
starburst model is consistent with all measurements for a 
bubble radius of r w = 1 h~ 1 Mpc and a volume filling factor 
of 2%. A bubble size of r w = 1 h _1 Mpc appears, however, 
still difficult to achieve for galactic winds. If for some reason 
the velocity shear around haloes is smaller than suggested 
by our (DM only) simulation this would further reduce the 
required bubble size and filling factor. The correlation length 
for our "starburst" model is 2.5 h _1 Mpc, slightly lower than 
the observations. Identifying the LBGs with satellites in- 
falling on high mass haloes, as in Kolatt et al. (1999), would 
further improve the agreement with the observed clustering 
properties of LBGs. However, this would be at the expense 
of matching the observed small-scale flux F, since the veloc- 
ity shear increases with halo mass. Therefore we advocate 
the model in which starbursts occur in small mass haloes. In 
a ACDM cosmology of normalisation erg = 0.9, this model 
yields a reasonable match to the observed clustering of LBGs 
and transmitted flux F. The mean transmitted flux of our 
simulated spectra at the smallest separation is 30 percent 
smaller than that of the ASSP sample. We would thus ex- 
pect the observed mean transmission to decrease substan- 
tially for larger samples if the model is correct. A sample 
with 12 LBGs at separation s < 0.5 h~ Mpc should reduce 
the error to AF ~ 0.1. This should help to discriminate 
against alternative models like the suggestion of Croft et al. 
(2002) that LBGs are located preferentially in low-density 
environments. 
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APPENDIX A: EVOLUTION OF BUBBLES 

To explain the lack of Hi close to the haloes, the dilute, 
highly ionised bubble has to survive for a sufficiently long 
time. In this appendix we present general order of magnitude 
estimates of the wind/shock energetics, relate that to the 
observed speeds, and estimate cooling time scales. We heav- 
ily rely on the work of Sedov (1959), Weaver et al. (1977), 
Ostriker & McKee (1988) and Tegmark et al. (1993). 



Al Wind energetics 

The main characteristics of the wind, such as the tempera- 
ture of the inner plasma, the size r w of the bubble swept up 
by the shock and the wind velocity w w , depend on the total 
energy E w released in the outflow. Neglecting energy losses 
E„ can be estimated as 

1 



4.3 x 10 56 erg 



(Al) 



0.1 h^Mpc 



( \\ 

VlOOkms- 1 / 9 



where M w is the gas mass inside r w . To get the numerical es- 
timation we have assumed that M w = (47r/3)r w 3 /9(,A 9 where 
p b is the mean background density and A 9 is the density con- 
trast inside r w . The characteristic radius r w and velocity v w 
of the wind are poorly constrained, but recent observations 
suggest that r w JjO.l h _1 Mpc and w w JjlOO kms" 1 (Heck- 
man et al. 2000, Pettini et al. 2001, Adelberger et al. 2003). 
On scales r ~ 0.1 h _1 Mpc the value of A 9 is typically about 
2 — 5, which gives E w ~ 10 57 erg. This value is consistent 
with the amount of energy Esn released by supernovae dur- 
ing a burst of duration hurst ~ 10 7 yr and of star formation 
rate M* ~ 10 M Q yr _1 . 

To assess whether the gas will be expelled from the halo 
one can apply a simple binding argument (e.g. Dekel & Silk 
1986). On the one hand the potential energy U of the gas 
which lies within the halo is U oc MM s /vm, where tm is 
the characteristic radius of the halo, tm oc M 1 * 13 , and M s 
is the mass of the gas which lie within a radius tm- The 
mass of the gas in the halo is M g oc f+fbM, where is 
the star formation efficiency and /;, the average fraction of 
baryons in the halo. The potential energy is then U oc M 5 ^ 3 . 
On the other hand, if one assumes that winds are driven 
by supernovae, the energy released in the wind is E„ oc 
M+thurst, where M* oc M g is the star formation rate. Hence, 
according to this crude estimate, the ratio of the kinetic to 
binding energy is E w /U ~ M~ 2 ^ 3 . We therefore expect r w 
to be a function of M , and outflows to escape more easily 
from low mass haloes (e.g. from dwarf galaxies). On this 
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Figure Al. The cooling radius r coo ;(i) (solid curves) for various 
cooling time rate. The dashed curve is r s h(t). 



latter point, the measurements of Adelberger et al. (2003) 
seem to indicate that, if outflows are the cause of lack of 
Hi absorption in the observed transmitted flux F, they can 
escape from high mass galaxies as well, and affect the IGM 
properties on comoving scale as large as 0.5 h _1 Mpc. 



A2 Propagation of adiabatic shock 

A wind ejects energy in the IGM, driving a shock which 
heats up and collisionally ionises the material it encounters. 
We briefly review now the evolution of the shock front and 
the ionised bubble it encompasses, before radiative losses be- 
come important. We assume that the gas follows the perfect 
gas law P = (7 — l)up, where P, u and p are, respectively, 
the pressure, the internal energy and the density of the gas, 
and 7 is the ratio of specific heats, which is 7 = 5/3 for a 
mono-atomic gas. If the expansion of the universe can be 
ignored, the shock is well described by a strong adiabatic 
blast wave which propagates into the IGM . In this regime, 
accurate self-similar solutions exist. Let r s h(t) be the posi- 
tion of the shock front at time t after the burst. The jump 
conditions relate the density pigm of the unperturbed IGM 
to the density p s h, the pressure P a h and the temperature T s h 
at the shock front. For a strong adiabatic shock the jump 
conditions for pressure and density can be expressed as 



(7 + 1) „ 
Psh = ( 7 -l) pIGM ' 



P sh = -^-Pigm^ 2 . (A2) 
7 — 1 at 



Hence, p Bb — 4piGM for an adiabatic index 7 = 5/3. As- 
suming that the flow is self-similar, r s h can be expressed 
as a function of the similarity variable £. The position of 
the shock front corresponds then to a fixed value of £ ~ 1. 
Hence, r s h is given by (Sedov 1959) 
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where E„ ~ 10 57 erg is the energy released in the wind. 
Inserting the value of _E W in the expressions of r s h and T s h, 
we find 
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Since the temperature behind the shock front T oc P/p is 
a decreasing function of r, T reaches its minimum at the 
shock front. This solution holds only in the regime where the 
expansion is negligible, i.e. t < ijj where tn is the Hubble 
time scale. Self-similar solutions for t > tn including this 
latter complication exist (Ostriker & McKee 1988), but since 
energy and momentum conservation break down for t ~ tH , 
it is difficult to express the asymptotic solution as a function 
of the initial E w . 



A3 Cooling time-scales and survival of bubbles 

The evolution of such an expanding shock depends critically 
on cooling which will lower the temperature of the shell and 
its interior, thereby increasing the amount of Hi hydrogen 
available for Lya resonant scattering. The expanding shell 
cools due to radiative losses. The shock will undergo a phase 
of shell formation which will occur when the energy dissi- 
pated in a volume oc r s h(i) 3 will be of the order of E„. 

Cooling is expected to be more efficient at the shock 
front since, at r = r s h, the density (temperature) is larger 
(lower) than at r < r s h. Thus, shell formation first occurs at 
the shock front. As long as we are in the regime of low den- 
sity, p <C Pigm, and high temperature, T ^> 10 6 K, radiative 
cooling is dominated by thermal bremsstrahlung. However, 
for temperature T <^10 6 K, atomic line cooling might con- 
tribute significantly to the cooling rate A (in ergcm 3 s _1 ) 
if the wind is enriched with heavy elements. For a wind of 
metallicity Z ~ O.IZq (e.g. Adelberger et al. 2003), the cool- 
ing rate is A ~ 10~ 22 erg in the range 10 4 <T ^10 6 K (e.g. 
Sutherland & Dopita 1993). 

We assume that cooling at the shock front, where the 
temperature is T <10 6 K, is dominated by atomic line tran- 
sitions. This is true when the temperature of the shell is 
T s h <^10 6 K, i.e. for times t ;>10 8 yr. The cooling time scale 
at the shock front is then t coo i(r B h, t) — 3fcsT , s h/(2n sh A), 



tcooi (^'sh, £) — 



1.9 x 10 10 yr 



VL b h 2 
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10 7 yr 
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where the cooling rate A22 is in unit of 10 
and is assumed to be constant in the temperature range 
10 4 ^T sh ^10 6 K. 

Once the shell starts forming, the temperature T sh of 
the shell cools to ~ 10 4 K on a very short time scale. How- 
ever, in the range T <10 4 K, the cooling rate A drops sharply 
from ~ 10~ 22 down to ~ 10~ 26 erg cm 3 s _1 (Sutherland & 
Dopita 1993) and the shell cannot cool below temperature 
T ~ 10 4 K. We define now a radius r coo i(t) such that the 
cooling time scale t coo i as a function of radius r satisfies 
tcooi < t for r > r coo i. According to this definition, the 
gas between r coo l{t) and r s h(t) has cooled onto the shell 
by the time t. To express t coo i as a function of r and t, it 
should be noted that the number density n w (7", i) behind the 
shock front (in cm ) scales as n w (r, i) = (r/r s h) a n B h(t) 
with a ~ 4.5 typically. This relation is however valid in 
the asymptotic regime only. Therefore, assuming that the 
plasma has reached pressure equilibrium, the temperature 
T w (r,t) behind the shock front is T w = (r sh /r) Q T sh , and the 
radiative cooling time is 



t C ooi(r,t) = (-jr) t coo i(r Bh ,t) 



(A7) 



On Fig. Al the solid curves are r coo i(t) for various cooling 
time rate A22, and the dashed curve is r s h(t). The shell starts 
forming at the shock front when r coo i(t) intersects r sh (i). For 
larger time, cooling occurs for r coo i(t) < r < r s h(t) (shaded 
area). Since most of the gas swept up by the wind lies at 
r <>' s h, a large fraction of gas cools down to T ~ 10 4 K 
shortly after the shell formation. For A22 = 1, about 30% 
of the mass enclosed in the bubble has cooled into the shell 
over a Hubble time tn, where tn — 1.2 x 10 9 (f2 m ft 2 ) 
at z = 3. For a high metallicity wind Z JjO.IZq, cooling 
becomes severe for t J> a few 10 yr, a time scale shorter 
than the Hubble time scale tH at z = 3. On the contrary, 
for a low metallicity wind, Z <C O.IZq, the cooling rate is 
A22 <C 1 and cooling will not become significant for t <tH- 
Hence, most of the bubbles will not survive over a Hub- 
ble time scale unless the cooling time rate A22, or the wind 
metallicity Z, is very low. For radii r <r coo i(t), the tem- 
perature is higher than 10 6 K, and thermal bremsstrahlung 
and inverse Compton cooling prevail over line cooling, but 
the corresponding cooling times are still smaller than the 
Hubble time at z = 3. 



A4 Absorption by the cold shell 

To compute the thickness As(t) of the shell, we assume that 
the pressure is constant through the inner plasma and the 
shell. Since the temperature of the shell, T s h ~ 10 4 K, is 
much lower than the overall temperature T w J>10 6 K of the 
highly ionised plasma, we expect As(t) <C r s h(t), i.e. the 
shell is thin. Taking into account this assumption, we find 
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4,2 a result consistent with a thin shell approximation. Assum- 
ing T w = 100T s h = 10 6 K and neglecting the expansion of 
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the Universe gives As <0.01r sh . As should thus not exceed a 
few h _1 kpc, a scale which cannot be resolved in our simula- 
tion. The optical depth for resonant Lya scattering is given 
by 



<7on s hAs 
1.6t?A s 



(A9) 
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0.1 h^Mpc 



where r\ ~ 1 is the Hi neutral fraction at temperature 
T ~ 10 4 K. Lya absorption is expected for reasonable val- 
ues of the density contrast Ag ~ 2 — 3, the typical comoving 
radius of the bubble r w ~ 0.1 h _1 Mpc, and the baryon frac- 
tion Q, b h 2 ~ 0.02. However, since we cannot resolve the thin 
shell in our simulation, we have nevertheless ignored its Lya 
absorption. 



APPENDIX B: NUMERICAL CONVERGENCE 

To assess the sensitivity of the transmitted flux to numer- 
ical resolution, we perform a calculation with the highest 
resolution simulation of the series, which has a mass res- 
olution m — 1.66 x 10 8 M Q /h. The halo catalog and the 
flux in the simulation are computed as in Section §3 with 
a pixel resolution A = 60 h _1 kpc (comoving). We added 
a Gaussian deviate of variance a v = 150 kms -1 to the 
halo positions. In Fig. Bl we show the transmitted flux as 
computed in the "starburst" scenario for a bubble radius 
r w = 1 h _1 Mpc. The mass resolution is m = 9.52 (dashed 
curve) and m = 1.66 x 10 8 Mg/h (solid curve). In the lat- 
ter the strong absorption lines are better resolved, thereby 
slightly reducing the normalization of the optical depth A. 
Therefore, the optical depth of the absorbing material is 
lower, and the transmitted flux is slightly larger, by about 
AF ~ 0.02 on scales of 2 < s < 4 h _1 Mpc. 




Figure Bl. The transmitted flux F in the "starburst" scenario 
with bubbles of radius r w = 1 h Mpc, with DM particle mass 
resolution m = 9.52 X 10 8 (dashed curve) and 1.66 X 10 8 M /h 
(solid curve). The horizontal dashed line is the flux normalization, 
(F) = 0.67. 



